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ABSTRACT 



The in-line force acting on smooth circular cylinders 
placed in oscillatory flow has been measured using a recently 
constructed U-tunnel. 

The drag and inertia coefficients have been determined 
through the use of the Fourier analysis and found to depend 
on Reynolds number as well as Keulegan-Carpenter number for 
Reynolds numbers about greater than 25,000. A new parameter 
namely the "Frequency Parameter" has been defined as the 
ratio of the Reynolds number to the Keulegan-Carpenter number, 
and the existence of unique relations between the frequency 
parameter and the drag and inertia coefficients have been 
shown . 

It is recommended that the experiments be extended to 
even larger Reynolds numbers using the same experimental 
technique and apparatus. 
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NOMENCLATURE 



A Virtual amplitude of the motion 

A^ Amplitude of the motion 

a Maximum acceleration 

m 

C D Average drag coefficient 

Average inertia coefficient 
Maximum in-line force coefficient 

FMAX 

D Diameter of the test cylinder 

F Instantaneous total force acting on the test cylinder 

F d Drag force acting on the test cylinder 

Fj Inertia force acting on the test cylinder 

F Measured force 

m 

g Gravitational acceleration 

H Distance between pressure taps and mean water level 

K Keulegan-Carpenter number 

L Length of the test cylinder 

Re Reynolds number 

s Distance between the pressure taps 

T Period of oscillations 

t Time 

U Maximum velocity 

m J 

u Instantaneous velocity 

w Width of the test section 

£ Frequency parameter 

y Specific weight of fluid 
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A Differential pressure 

0 Phase angle 

v Fluid kinematic viscosity 

p Fluid density 
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I. 



INTRODUCTION 



Information about the time-dependent forces acting on 
bluff bodies in general and circular cylinders in particular 
has considerable practical interest in ocean and wind engi- 
neering, as well as in basic understanding of fluid mechanics. 
Extensive discussion about the flow-induced forces and oscil- 
lations exists in the literature, but the basic hydrodynamic 
data are lacking, particularly at higher Reynolds numbers 
which are of current practical interest. 

Much of the present knowledge on the hydrodynamics of the 
oscillatory flow about the cylinders has been obtained by 
means of model tests in wave channels at Reynolds numbers 
generally two to three orders of magnitude smaller than proto- 
type Reynolds numbers. These model tests, which have been 
primarily concerned with in-line forces, have disclosed that 
the drag and inertia coefficients may be correlated with the 
amplitude of the motion relative to a characteristic dimension 
of the body and that the effect of the Reynolds number is 
obscured by that correlation. The field data have shown 
rather large scatter and no systematic correlation with either 
the relative amplitude or the Reynolds number. Understandably, 
there is considerable controversy in the literature regarding 
the effect of viscosity, free surface, etc. 

In view of the foregoing considerations, the present' 
research program was undertaken with the following objectives: 
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(a) To furnish data, obtained under carefully controlled 
laboratory conditions, about the in-line and trans- 
verse forces acting on smooth cylinders in a 
harmonically oscillating fluid at Reynolds numbers 
up to 400,000 and 

(b) To identify the physical mechanism and parameters 
responsible for the correlation or scatter of the 
force-transfer coefficients. 

This thesis deals only with the in-line forces acting on 
smooth circular cylinders. It does not deal with ocean waves, 
non-harmonic fluid oscillations, wave and current combinations 
and its consequences, diffraction effects, free-surface and/or 
wall proximity effects, fluid elasticity or hydro-elasticity 
of flexible or flexibly supported cylinders in harmonic motion, 
or finally the interference effects between neighboring 
structural elements. Furthermore, no attempt is made to offer 
a chronological and/or critical survey of the fluid loading 
on cylinders or offshore structures. Fairly complete accounts 
in the context of wave forces are given by Wiegel (1) , 

Hogban (2) , and Grace (3) , where an extensive list of refer- 
ences can be found. Only those works which have a direct 
bearing on the evolution and discussion of the present data 
will be reviewed, in some detail, wherever appropriate. 

In closing the introduction it should be emphasized that 
the force-transfer coefficients must be obtained experimentally 
for the analysis of the harmonic flow about cylinders may 
never be forthcoming. 
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II. ANALYSIS OF THE DATA AND THE FORCE COEFFICIENTS 

A completely satisfactory analysis of the resistance in 
unsteady separated flow has escaped the concentrated efforts 
of many researchers. No theoretical model can, at present, 
predict the complete force and flow characteristics of a 
harmonic motion about a circular cylinder. In the absence of 
such an analysis, the most serious difficulty lies in the 
description of the time-dependent force. Other difficulties 
arise in the description and interpretation of the history 
of the flow and the effect of the vortices and turbulence in 
the flow approaching the body at a particular time in the 
oscillation. One approximate and physically meaningful way 
around these difficulties is to assume, following Stokes' 
classic analysis of the oscillating pendulum, that the total 
time-dependent in-line force may be expressed as a sum of the 
velocity-squared dependent drag and the acceleration dependent 
inertial force, each with a suitably averaged force coefficient. 
This then is the basis of the so-called Morison's (4) equation. 
For a circular cylinder it may be written as 

F = Fp + Fj = 0.5 C d LD p|u|U + 0.25 tt p L D C M ^ (1) 

For an oscillating flow represented by U = - U^CosG , with 
0 = 2rrt/T , the Fourier averages of C D and are given 

by Keulegan and Carpenter (5) as 
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in which F^ represents the measured force, L and D the 
length and the diameter of the cylinder, U m the maximum 
velocity and, T the period of the flow oscillations. Evi- 
dently, C D and represent the first two terms in a 

complete series expansion of the normalized force in terms of 
the odd integers in sines and cosines. Previous studies by 
Keulegan and Carpenter (5) and Sarpkaya and Tuter (6) have 
shown that C D and as given by equations (2) and (3) 

are the most significant ones and sufficient to represent the 
measured force adequately. 

A. GOVERNING PARAMETERS 

The coefficients cited above will have to be correlated 
through the use of suitable parameters in order to show that 
they have some degree of universality within the range of the 
parameters encountered. A simple dimensional analysis of the 
flow under consideration (uniform harmonic motion about a 
circular cylinder placed with axis normal to the flow) shows 
that the time dependent force coefficient may be written as: 

= f (U T/D , U D/v , t/T) • (4) 

0.5 p U D L 
m 

= f (K, Re, t/T) (5) 
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Evidently, U^T/D may be replaced by 2irA/D or simply by 
A/D where A represents the amplitude of the oscillations. 

There is no simple way to deal with equation (5) even for 
the most manageable time-dependent flows. The evaluation of 
the instantaneous values of the force coefficients is not at 
present feasible. Another and perhaps the only alternative 
is to eliminate time as an independent variable in equation 
(5) and consider suitable time-invariant averages of the 
force coefficients. Thus, one has 



'D 



'M 



'FMAX 



= f. (U T/D, Re) 
1 m ' 



( 6 ) 



Even this equation (6), as simple and as idealized as it is, 
gives rise to many questions: Do the averaged force co- 

efficients really depend on both K and Re ; can one obtain 
meaningful conclusions by plotting the data for a given co- 
efficient with respect to, say, K and connecting points 
having equal Re or vice versa ; how should the experiments 
be conducted so that equation (6) yields manageable plots ; 
which of the two parameters has a more pronounced effect on 
the force coefficient under consideration ; why has there 
been considerable scatter (see, e.g., Wiegel (1)) in the field 
data when plotted with respect to either K or Re ; are 
there ranges of K and Re in which the effect of one is 
obscured by a reasonable correlation of the force coefficients 
with the other? These and similar questions have been raised 
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by many other investigators and attempts were made to establish 
suitable correlations. The state of the art is such that the 
past conclusions and conjectures can be critically scrutinized 
only through the acquisition of data of high intrinsic quality 
obtained under controlled laborabory conditions with relatively 
simple harmonic flow situations. The purpose of such an effort 
is by no means to remove the need for actual full scale 
experiments. In fact it is to encourage full scale experiments 
and to enable those concerned to interpret and better understand 
the factors effecting the force-transfer functions. 

Past experience (5,6) has shown that the force coefficients 
obtained under controlled laboratory conditions are primarily 
functions of K at relatively small Reynolds numbers 
(Re < 25,000) and that the effect of viscosity is obscured by 
the fairly good correlation between K and the force co- 
efficients. Again, previous efforts and the reasoning based 
on dimensional analysis have shown that there is an undeniable 
effect of the Reynolds number. Thus, means have to be devised 
to delineate the effect of both K and Re or some other 
viscosity dependent parameter. 

It is apparent that the maximum velocity appears in 

both K = U T/D and Re = U D/v . Simple rules of dimensional 
m m 

analysis state that one obtains the maximum amount of experi- 
mental control over the dimensionless variables, if the 
original variables, that can be regulated, each occur in only 
one dimensionless product. Thus, if U m is easily varied 
experimentally, then U m should occur in only one of the 
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independent dimensionless parameters. With this hint in mind, 

2 

Re in equation (6) is replaced by Re/K = D /vT . This para- 
meter shall be called the "Frequency Parameter" and denoted by 
6 , so that: 



Evidently, for a series of experiments conducted with a 
cylinder of a given diameter D in water (of uniform and 
constant temperature and density) undergoing harmonic oscil- 
lations with a constant period of T , 3 is held constant, 

then the variation of a force coefficient with K may be 
plotted for constant values of & . Subsequently, one can 
easily recover the Reynolds number from 

Re = BK (8) 

and connect the points , on each 8 = constant curve repre- 
senting a given Reynolds number for a family of suitably 
selected values of the Reynolds number. Such a procedure 
eliminates the difficulty of trying to draw contours of con- 
stant K , or constant Re , or constant C D or versus 

K or Re , or K versus Re . Suffice it to note that the 
data reported here-in shall be analyzed according to the 
relationship 



C. 

l 

a coefficient = f^ (K,B) • (9) 
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and the Reynolds number will be used in the manner cited above. 
The power of this new procedure (new, as far as the wave force 
analysis is concerned) will become apparent later. 
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III. EXPERIMENTAL EQUIPMENT AND PROCEDURES 



A. U-SHAPED OSCILLATING-FLOW TUNNEL 

Experiments carried out by Sarpkaya and Tuter [6] in the 
initial phases of the study with small smooth cylinders at 
relatively low Reynolds numbers have proved the versatility 
and the usefulness of a U-shaped oscillating-f low apparatus. 
Thus , in an attempt to achieve larger Reynolds numbers , it 
was only natural to construct a larger U-shaped tunnel. 

Among the various designs made by Professor Sarpkaya, the 
one shown in figure 1 was finally selected for construction. 

A photograph of the completed and fully instrumented structure 
is shown in figure 2. It consists of nine modules for ease 
of construction, transportation, and mounting. Each module 
is made of 0.95 cm. (3/8 inch) aluminum plates and reinforced 
1.27 x 10 x 46 cm. (1/2 x 4 x 18 inches) aluminum flanges 
welded to the plates. The modules were assembled with the 
help of an air drying silicon rubber seal between the flanges 
of two adjacent modules and 2.54 cm. (1 inch) steel bolts 
placed 15 cm. (6 inches) apart. The inside of each module 
was precision machined so that the largest misalignment was 
about 1 mm. (0.04 inches). 

Prior to the description of its instrumentation and 
operation, a few words are necessary about the general shape 
of the tunnel. The cross-section of the two legs is 
183 x 91.5 cm. (6x3 feet). This selection was dictated 
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Figure 1. Schematic Drawing of the U-Tunnel 




Figure 2 . U-Tunnel 
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Figure 3. Schematic Drawing of the 
Butterfly Valve System 
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by several considerations such as available ceiling height, 
pressures to be encountered and hence the structural and 
economic considerations, desire to obtain a virtual amplitude* 
or velocity of oscillation at least twice that of the free 
surface, period of oscillation, Reynolds number and the 
relative amplitude A/D desired, natural-damping of the 
oscillations, and the magnitude and the frequency of the 
forces. The length of the horizontal test section was chosen 
larger than twice the virtual amplitude to insure fully de- 
veloped uniform flow at the test section. Finally the two 
corners of the tunnel were carefully streamlined to prevent 
flow separation (see figure 1) . This design proved to be 
more than adequate for no separation was encountered, and 
also the desired frequency and amplitude of oscillation were 
achieved. 

The auxiliary components of the tunnel consisted of 
plumbing for the filling and emptying of the tunnel, butterfly- 
valve system, and the air supply system. The plumbing 
consisted of simple piping, valves, a small pump, and a 
filter. None of these will be described further. 

The butterfly-valve system (mounted on top of one of the 
legs of the tunnel) consisted of four plates, each 4.6 cm. 

(18 inches) wide and 92 cm. (36 inches) long. A 2.54 cm. 

*The virtual amplitude is defined as the amplitude of oscil- 
lation which the cylinder-experiences. Here it is exactly 
twice the amplitude of oscillation of the free surface in 
one leg of the tunnel. 
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(1 inch) steel shaft was placed at the axis of each valve 
plate. (see figure 3) Aluminum housings supported both ends 
of the shaft with self-aligning ball bearings. A 15 cm. 

(6 inches) gear was attached to one end of each shaft which 
extended beyond the bearing. All four valve plates were then 
aligned and driven by a simple rack and pinion system. The 
rack was actuated by an air-driven piston with the help of 
a three-way valve connected to the laboratory air supply system. 

The valves, in their closed position, completely sealed 
the top of one of the legs of the tunnel. The top of the 
other leg was left open. Initially, the butterfly valves 
were closed and air was admitted to that side of the tunnel 
to create the desired differential water level between the 
two legs of the tunnel. Then the valves were opened quickly 
with the help of a pneumatically-driven three-way control valve. 
This action set the fluid in the tunnel in oscillatory motion 
with a natural period of T = 5.272 seconds. A series of 
experiments was conducted with one of the test cylinders to 
evaluate some of the experimental characteristics of the 
tunnel. It was found that the damping of the motion is such 
that the amplitude of oscillations decreases about 3 mm. 

(0.12 inches) per cycle and about 0.15 mm. (0.06 inches) per 
cycle for amplitudes half or less than of the maximum. Thus, 
over a period of four or five complete cycles of oscillation 
at any amplitude, the amplitude, velocity, and the ac- 
celeration of the fluid changed about one percent. Evidently 
the forced oscillations of the fluid, if such a method were to 
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be employed, cannot yield the amplitude to an accuracy better 
than one percent. Additionally in such a method one has to 
contend with some high frequency vibrations, however small, 
superimposed on the acceleration. These -result from the 
cyclic operation of the butterfly valves. It is because of 
these considerations that the experiments were carried out by 
letting the system damp out the amplitude over many cycles of 
oscillations. The advantages of the method adopted become 
apparent very quickly. Firstly, the oscillations were so 
smooth and quiet that one could not know or even hear that 
5000 gallons of water were in oscillation. The elevation, 
acceleration, and all force traces were absolutely free from 
secondary oscillations so that no filters whatsoever were 
used between the transducer output and the recording equipment. 
Secondly, the method adopted enables one to cover all possible 
values of K for a given 3 and see the evolution of the 
forces over a period of about 30 minutes. 

B. CIRCULAR CYLINDER MODELS 

Seven circular cylinders with diameters ranging in size 
from 20.3 cm. (8 inches) to 5 cm. (2 inches) have been used 
in this study. The cylinders were turned on a lathe from 
aluminum pipes or plexiglass rods and polished to a mirror- 
shine surface. There is no doubt that the surfaces were 
hydrodynamically smooth. The length of each cylinder was 
such that it allowed 0.08 cm. (1/32 inch) gap between the 
tunnel wall and each end of the cylinder. As will be noted 
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later, the cylinder was prevented from moving towards one or 
the other wall by means of small O-rings attached to the round 
cantilever end of the force transducers. A double-ball pre- 
cision bearing (SKF-2303-J) with an approximately 1.5 cm. 

(0.6 inches) bore was inserted at each end of the cylinder 
in aluminum housings which sealed the cylinder air tight. 

The other face of each bearing was flush with the end of the 
cylinder . 

C. FORCE MEASUREMENTS 

Two identical force transducers, one at each end of the 
cylinder, were used to measure the instantaneous in-line and 
transverse forces. The basic transducer was manufactured by 
B.L.H. Electronics, Inc. under the trade name LBP-1 and 
catalogue no. 420271. A photograph of the gage assembly is 
shown in figure 4. The gage had a capacity of 2224 N 
(500 lbf.) with an overload capacity of 200 percent. The de- 
flection of the gage under 500 lbf. load was 0.25 mm. (0.01 
inch) . For the largest cylinder and amplitude encountered 
in the experiments, the maximum load was less than 200 lbf., 
and the deflection of the beam was less than 0.2 mm. 

(0.008 inches ) . 

A special housing was built for each gage so that it can 
be mounted on the tunnel window and rotated to measure either 
the in-line or the transverse force alone. The bellows which 
protected the strain gages had to be waterproofed in such a 
manner that they would not adversely effect the operation 
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Figure 4. Gage Assembly 







of the gages when subjected to about 6 m. (20 feet) water 
pressure at temperatures 18 C (64 F) to 74 C (165 F) . For 
this purpose the bellows were completely filled with Dow 
Corning 3140 (RTV coating without bringing the rubber into 
contact with air during the filling operation. After filling, 
the ends of the bellows were sealed air tight with special 
clamps. The silicon rubber remained in its original liquid 
form throughout the operation of the gages. 

The cylinders were placed in the test section by re- 
tracting the gages from their housing and then pushing them 
into the bearings mounted on each end of the cylinders. As 
noted earlier the O-rings placed on the cantilever end of each 
gage prevented the test cylinders from moving sideways to- 
wards one or the other wall and helped to exactly set the 
0.8 mm. (1/32 inch) space between the cylinders and the 
tunnel wall. The cylinders were free to rotate at the 
application of a slight torque by hand. 

After the mounting of the first cylinder, the exact 
angular position of the gages within their housing had to 
be determined and set with a pin so that the gages measure 
either only the in-line or the transverse force. For this 
purpose, an approximately 400 N load was hung on the cylinder 
with a lubricated nylon rope. The in-line force (acting in 
the horizontal direction) was observed on the amplifier re- 
corder system. Then the gage was rotated in small increments 
until the in-line force was exactly zero. A final check was 
made by measuring the outputs of the gages with a precision 
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voltmeter. The position of each gage was marked and set with 
a pin. Finally, four bolts were placed on the gage housing 
to hold the gages rigidly in position. Removal of these 
bolts and the pin allowed the rotation of the gages exactly 
90 degrees, after which the bolts and pin were placed in 
position. In this manner the gages were capable of measuring 
either the in-line or the transverse force without any "cross 
talk" between the in-line and transverse forces. Ordinarily, 
one gage was set to measure in-line force and the other gage 
the transverse force. At times both gages were used to measure 
only the in-line or the transverse force. 

The calibration of each gage was accomplished by hanging 
loads in the middle of the cylinder after setting both gages 
to sense only the transverse (here vertical) force. The di- 
rectional sensititivity of the gages was also checked by 
applying identical loads upwards on the test cylinders with 
the help of a hook-cantilever arm attached to the top of the 
tunnel outside the test section. Repeated calibrations have 
shown that the gages are perfectly linear up to 2000 N; they 
yeild the same signal for loads applied either downward or 
upwards; and that the gages, together with the electronic 
system to which they were eventually connected were capable 
of sensing forces as small as 0.1 N (0.02 lbf.). 

The analog data appearing on a two channel Honeywell 
recorder in the form of either in-line force versus amplitude 
or transverse force versus amplitude were were read at every 
0.1 seconds. The evaluation of the data is discussed later. 
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D. ACCELERATION, ELEVATION, OR VELOCITY MEASUREMENTS 

It is because of the extreme importance of the accurate 
measurement of the instantaneous value of these quantities 
that they are discussed here separately. 

Firstly, it should be noted that the measurement of the 
amplitude, acceleration, elevation, or the velocity is a 
matter of interpretation of the signal received from the 
appropriate transducer in light of one of the following 
expressions 



n - 2 ttA _ dU 

m T ' a m dt 



il) ^ ^ _ 2 tt 
T ' A T m 



in which T = 5.272 seconds for the experiments reported here. 

Three transducers were used to generate three independent 
d.c. signals, each proportional to the instantaneous value 
of one of the quantities cited above. The first one con- 
sisted of a seven foot long platinum wire stretched vertically 
in one leg of the tunnel. The output of the capacitance-wire 
bridge was connected to an eight channel amplifier-recorder 
system. The response of the wire was found to be perfectly 
linear within the range of oscillations encountered. The wire 
was capable of yielding a measurable signal for changes in 
water elevation as small as 0.8 mm. (1/32 inches). Such a 
sensitivity was not, however, always desirable for the in- 
stabilities on the water surface gave rise to small oscillations 
in the analog records. The effect of such instabilities were 
practically eliminated by placing the wire along the axis of 
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30 cm. (1 foot) diameter and 213 cm. (7 feet) long thin 
plastic pipe. Be that as it may, the use of this wire was 
rendered unnecessary due to the use of a more reliable method. 

The second method consisted of the measurement of the 
instantaneous acceleration by means of a differential-pressure 
transducer connected to two pressure taps placed horizontally 
61 cm. (2 feet) apart and 91.5 cm. (3 feet) to one side of 
the test section. The output of the transducer was again 
connected to the eight channel recorder. The instantaneous 
acceleration was then calculated from Ap = ps(du/dt) where 
Ap is the differential pressure, s the distance between the 
pressure taps and dU/dt is the instantaneous acceleration 
of the fluid. The effect of the pressure drop due to the 
viscous forces over the distance s was calculated to be 
negligible . 

The third method again consisted of the measurement of 
the differential pressure between two pressure taps. The two 
taps were placed symmetrically on the two vertical legs of 
the tunnel at an elevation 96.52 cm. (38 inches), (see 
figure 5) below the mean water level, i.e., H = 38 in. 
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Applying Bernoulli's equation for unsteady flow between each 
pressure tap and the instantaneous level of water, it is easy 
to show that twice the amplitude of the free surface oscillation 
(virtual amplitude) is given by 



A = 2A = A P/l— , (10) 

l-i(^L) H 

g v T ' 

in which g and T are constant and H is kept constant. 

Thus the signal of this transducer yielded the virtual amplitude 
or the maximum velocity in each cycle. It was entirely free 
from noise or small free surface effect. The transducer was 
calibrated and its linearity checked before each series of 
experiments . 

Suffice it to say that all three methods gave nearly 
identical results and yielded the amplitude, velocity, or 
acceleration to an accuracy of about two percent relative to 
each other. These comparisons as well as the perfectly 
sinusoidal and noise-free character of all pressure and force 
traces speak for the suitability of the unique test facility 
used in this study. 

E. DATA REDUCTION 

Experiments were repeated at least three times for each 
cylinder. Often only two, but at times as many as four runs 
were evaluated. Suitably selected cycles were read every 0.1 
second and the tabulated data were transferred to the computer 
data cards together with additional relevant information such 
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as calibration factors, cylinder diameters, number of data 
cards, etc. A computer program calculated the force transfer 
coefficients previously discussed (see Appendix A) . 

Sample traces of the in-line force and the elevation are 
shown in figure 6. 
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Figure 6 . Elevation and Drag Force Traces 
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IV. DISCUSSION OF RESULTS 



A. BLOCKAGE AND LENGTH- TO-D I AMETER EFFECTS 

Attempts to achieve as high Reynolds numbers as possible 
in conducting wind-tunnel or water-tunnel experiments in- 
variably give rise to wall-interference effects which, of 
course, influence whatever measurements are desired. There 
are several blockage correction formulas for steady flows 
which might be used so that the wall-interference effects 
might be minimized. Unfortunately none of these formulas 
could be used in the present study for no one has demon- 
strated that the blockage effects in oscillatory flow are 
identical to those experienced in steady flows. 

The blockage ratio D/w , where D is the diameter of 
the cylinder and w the width or height of the test section, 
and the length- or span-to-diameter ratio, L/D , for the 
cylinders used in the present study are tabulated below 



II 

kO 


44 cm 


. (3 feet) , L = 

D 


90.885 

D/w 


cm . (2.9818 

L/D 


feet) ) . 
D 2 /vT 


20.244 


cm. 


(7.970 


inches) 


0.22 


4.49 


7968.7 


16.447 


cm. 


(6.475 


inches ) 


0.18 


5.52 


5259.9 


15.177 


cm. 


(5.975 


inches) 


0.17 


5.99 


4480.2 


12.674 


cm. 


(4.990 


inches) 


0.14 


7.17 


3123.2 


10.103 


cm. 


(3.978 


inches ) 


0.11 


8.99 


1985.2 


7.544 


cm. 


(2.970 


inches) 


0.082 


12.05 


1106.6 


6.349 


cm. 


(2.4996 


i inches) 


0.069 


14.31 


783.8 


5.057 


cm. 


(1.991 


inches ) 


0.055 


17.97 


497.2 


c 

II 

• 


05 x 


10 5 ft 2 /sec (0 


.009755 


stokes ) 
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Also shown in the above table are the $ = D /vT values for 
later use. Achenbach (7) and in some of the experiments of 
Fage and Falkner (8) the blockage ratios were 0.166 and 0.185 
respectively. The length-to-diameter ratio in Fage and 
Warsap's (9) experiments was 20.2 or 7.88, depending on the 
diameter of the two cylinders they used, as compared to 3.33 
in the experiments of Achenbach. It is generally observed 
that values of C D are smaller for cylinders with larger 
length-to-diameter ratio. Thus the presence of both the gaps 
and the larger L/D ratios could result in lower drag 
coefficients. Furthermore, the wake of the cylinder, near 
the cylinder ends is supplied with high pressure fluid from 
the front and as a result smaller values of C D are expected 
since the base pressure is increased over the value it would 
otherwise attain. In the subcritical range of Re, however, 
these effects appear to be negligible. In fact, Mosbach (10) 
found that in the subcritical range there is no effect of 
length-to-diameter ratio. 

In the present experiments, the gap cannot be eliminated 
by extending the cylinder into a cylindrical cavity within the 
two windows supporting the gages and the cylinder because of 
the fact that eight cylinders of different diameters were used. 
Obviously, it would have been too costly to build a window 
for each cylinder. It is believed that the very small gap 
used plus the cantilever end of the gage extending into the 
bearing in the cylinder minimized the supply of high pressure 
fluid into the wake from the front of the cylinder during one 
half of the cycle. 
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Returning to the discussion of the blockage effect it 
must be emphasized that the formulas used for steady flow 
correction effects cannot be applied to oscillating flows 
and that there is not a unique blockage correction for the 
entire period of the harmonic flow. This is evident from 
the fact that within a given cycle the fluid undergoes varying 
accelerations and velocities as the wake width, momentum 
deficiency, and the wake pressure change accordingly. Thus, 
a blockage correction made for the instant of maximum 
velocity is not applicable to the instant at which the maximum 
acceleration occurs. In view of the fact that there are no 
previous investigations, a series of experiments had to be 
conducted to determine the role of blockage in the flow under 
consideration. For this purpose a differential pressure 
transducer was connected to two pressure taps on the same 
side of the tunnel wall. One of the taps was placed on the 
wall directly above the axis of the test cylinder. The other 
tap was placed 76 cm. (29.92 inches) to one side of the first 
tap along a line parallel to the flow. A series of experiments 
were carried out with the 16.447 cm. (6.475 inches) cylinder 
and the differential pressure was recorded and compared with 
the differential pressure obtained from the acceleration 
transducer. Furthermore, to simplify the comparison both 
transducers were calibrated so as to render exactly the same 
output under identical calibration loads. The results have 
shown that the two differential pressures were nearly identical 
and that they were certainly within three percent of each other. 
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Often the two traces of two transducers were undistinguishably 
coincident. This somewhat surprising result is a clear 
indication of the fact that the blockage effect in harmonic 
flows is negligible at least for D/w ratios less than 0.18. 
Although no special attempt was made to interpret the lack 
of blockage effect in such flows it is believed that the 
presence of vortices on both sides of the cylinder together 
with the high periods of acceleration and velocity render 
the flow relatively more uniform at short distances away from 
the cylinder in the test section. Therefore, for the reasons 
cited above no blockage-effect corrections were applied to 
the data presented here. It might be of interest to note 
that had the flow been assumed uniform and had the maximum 
velocity for the largest cylinder and the Reynolds number 
were used to calculate a blockage effect correction through 
the use of one of the existing formulas, one would have 
found that such a correction would have amounted to about 
six percent. 

B. IN-LINE FORCE DATA 

In-line force coefficients calculated as previously de- 
scribed are plotted as a function of the so-called Keulegan- 
Carpenter number denoted here by K . A representative plot 
of C D versus K for the 16.447 cm. (6.475 inches) cylinder 
is shown in figure 7. It is evident from this figure that 
there is practically no scatter in the data even though it 
represents the result of four independent runs. This is 
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once again indication of the fact that the tunnel and its 
instrumentation are extremely reliable and that the data when 
evaluated with extreme regard relative to the phase angle are 
perfectly repeatable. As noted earlier the Reynolds number 
along this plot varies with K and $ . One can designate 
special values of the Reynolds number along this curve simply 
by writing Re = K£ or by finding the K values which 
correspond to a particular value of Re by writing K = Re/£ 
Prior to doing so mean lines have been drawn through the data 
as carefully as possible for each cylinder data. Then round 
numbers such as 30,000, 40,000, etc. were chosen for Re and 
the corresponding K values were calculated and indicated on 
each plot. This simple procedure resulted in figure 8 in 
which Cp is plotted as a function of K for all cylinders 
and the identical Reynolds number points are suitably 
connected. This procedure has also been repeated for the 
inertia coefficient and the result is presented in 

figure 9. It is evident from figures 8 and 9 that there is 
a remarkable correlation between the force coefficients, 
Keulegan-Carpenter number and the Reynolds number. The 
smoothness of the constant Reynolds number lines is another 
indication of the consistency of the data from one cylinder 
to another. 

Figure 8 shows that for a given Reynolds number, for 
Reynolds numbers larger than approximately 50,000, 
increases at first with increasing K , reaches a maximum, 
then decreases and then continues to rise gradually. For 
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Figure 7. Drag Coefficient Versus Keulegan-Carpenter Number 
for D = 6.475 Inches Cylinder 




40 



Figure 8. Drag Coefficient Versus Keulegan-Carpenter Number 
for All Test Cylinders 



Reynolds numbers larger than about 50,000 the initial rise 
in C D with K is very small or non existent. For very 
large Reynolds numbers and K , C D attains a practically 
constant value of about 0.5. For a given K , however, 
always decreases with increasing Reynolds number. It is 
evident from figure 8 that the change in with Reynolds 

number for Re smaller than approximately 25,000 is almost 
negligible. It is for this reason that Keulegan and Carpenter 
(5) and Sarpkaya and Tuter (6) reached the conclusion that 
C D does not depend on Reynolds number. The variations of 
C D with Re become most apparent between Re = 25,000 and 
Re = 150,000. It is also noted that for large values of 8 
such as B = 5260 and 4480 the values do not significantly 

change with Re . In other words C D values corresponding 
to B > 5000 yield the lower limit of the C D values whereas 
C D values corresponding to 8 - 500 consitute the upper 
limit of the C D values. The most significant changes in 
C D with K and Re occur between these two B values. 

Figure 9 shows the variation of the inertia coefficient 
with K and Re. The most significant changes in 
occur in the neighborhood of K - 12 . It is also around 
this value of K that changes significantly with Re . 

As discussed in connection with the variation of the drag 
coefficient, B values of 500 and 5000 more or less de- 
termine the lower and the upper limits of the inertia 
coefficient. Figure 9 also shows that for a given K 

increases with increasing Reynolds number and varies very 
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Figure 9. Inertia Coefficient Versus Keulegan-Carpenter Number 
for All Test Cylinders 



little for Re smaller than approximately 25,000. It is 
for this reason that in the studies by Keulegan and Carpenter 
(5) and Sarpkaya and Tuter (6) was thought to be inde- 

pendent of Reynolds number. The present results show that 
this is not so and that C„. varies with both Re and K 
throughout the K and Re range encountered in the present 
investigation. It appears that for very large Reynolds 
numbers C will approach a value of about 1.85. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The present investigation of the in-line forces acting on 
smooth cylinders in the range of K values from about 3 to 
200 and Re from 10,000 to 400,000 warranted the following 
conclusions . 

1. The drag and inertia coefficients could be determined 
quite accurately through the use of the present 
experimental system. 

2. Both the drag and inertia coefficients depend on the 
Reynolds and Keulegan-Carpenter numbers. For a given 
K , the drag coefficient decreases with increasing 
Reynolds number whereas the inertia coefficient 
increases with increasing Reynolds number. 

3. The most significant changes in both C^ and C 

occur in the neighborhood of K = 12 . 3 values of 

about 500 and 5000 determine the upper and lower 
limits of the two coefficients. 

4. Neither C D nor C M vary appreciably with Re for 
Re < 25,000. 

5. It is recommended that the experiments be carried out 
at even larger Reynolds numbers through the use of 
the existing systems by decreasing the viscosity of 
the water through heating and with cylinders of 
various artificial roughnesses. 
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6. It is also recommended that additional investigations 
such as the hydroelastic oscillations of flexibly 
mounted cylinders in harmonic flow and the hydro- 
elastic oscillations of ridigly mounted flexible 
cylinders be carried out. Finally the determination 
of the effect of wall proximity on the in-line and 
transverse forces acting on cylinders and other bluff 
bodies will be of extreme importance. 
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